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We examine the probabihty distribution function (pdf) of energy injection rate (power) in nu- 
merical simulations of stationary two-dimensional (2D) turbulence in the Lagrangian frame. The 
simulation is designed to mimic an electromagnetically driven fluid layer, a well-documented system 
for generating two-dimensional turbulence in the laboratory. In our simulations, the forcing and 
velocity fields are close to Gaussian. On the other hand, the measured PDF of injected power is 
very sharply peaked at zero, suggestive of a singularity there, with tails which are exponential but 
asymmetric. Large positive fluctuations are more probable than large negative fluctuations. It is 
this asymmetry of the tails, which leads to a net positive mean value for the energy input despite 
the most probable value being zero. The main features of the power distribution are well described 
by Craig's XY distribution for the PDF of the product of two correlated normal variables. We show 
that the power distribution should exhibit a logarithmic singularity at zero and decay exponentially 
for large absolute values of the power. We calculate the asymptotic behavior and express the asym- 
metry of the tails in terms of the correlation coefficient of the force and velocity. We compare the 
measured pdfs with the theoretical calculations and briefly discuss how the power pdf might change 
with other forcing mechanisms. 

PACS numbers: 47.27.Gs 



INTRODUCTION 

Since turbulence is an intrinsically dissipative phe- 
nomenon, it requires an external source of energy to sus- 
tain it. For a turbulent flow in a statistically stationary 
state, the rate of injection of energy into the system from 
this external energy source, the input power, is equal, on 
average, to the rate of dissipation of energy by small scale 
viscous processes. The fact that equality holds only on 
average is crucial. Since the input power is typically cal- 
culated as a product of an external force with the fluid 
velocity, both being fluctuating quantities, it is itself a 
fluctuating quantity with a full statistics of its own. Lo- 
cally in space or in time this injected power need not 
balance the corresponding rate of dissipation. In fact it 
need not even remain positive. 

Understanding the rate of energy injection in turbu- 
lence is of considerable importance in an engineering con- 
text, since it relates to the power required to overcome 
turbulent drag to sustain the rotation of a fan or turbine 
at a given speed in a turbulent flow (see, for example, [l|). 
In this context, much research focusses on the mean value 
of the power and how it scales with the Reynolds number. 
Such measurements, focussing on the mean power, in the 
context of Taylor-Couette flow were first performed by 
Lathrop et al. Subsequent research has focussed 

on measuring the pdf of the power fluctuations as well 
as the mean value in several turbulent systems includ- 
ing von-Karman flows [3, 0] , electro-convection [ll , wave 
turbulence ^ and turbulent convection . In this latter 



case, it was the heat transfer rather than the power which 
was measured. Much of the recent work has focused on 
finding macroscopic non-equilibrium systems on which 
to test various non-equilibrium "Fluctuation Relations" 
0, S [3] which have attracted much theoretical inter- 
est in recent years. In simple terms, these relations ex- 
press a symmetry of the probability distribution of some 
time-integrated quantity associated with the dissipation 
or entropy production in a non-equilibrium system. For 
a review, see [ll| and the references therein. 

All of the works cited above have focused on measure- 
ments of the statistics of global quantities. All have 
observed non-trivial probability distributions, typically 
with exponential tails but, to date, most of the discus- 
sions of the power distribution have been qualitative. In 
this article we use the notion of Lagrangian turbulence 
as a local measure of the energy injection into a tur- 
bulent flow. We measure the probability distribution 
from numerical simulations and give a quantitative ex- 
planation of the observed distribution in terms of Craig's 
XY-distribution for the product of correlated Gaussian 
variables. Although one can test the Fluctuation Rela- 
tion for Lagrangian power fluctuations this is not 
the purpose of the present article. We shall focus on the 
statistics of the power itself rather than it's time integral. 

The outline of the paper is as follows. We first give 
some background on turbulence in two dimensions, ex- 
plaining the difference between the direct and inverse 
cascades. We then give some details of our numerical 
simulations and introduce the notion of Lagrangian mea- 
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surements and the concept of Lagrangian power as a di- 
agnostic of the energy injection into a turbulent flow. We 
then derive some asymptotic properties of the probabil- 
ity distribution of the product of two correlated normal 
variables (Craig's XY-distribution). We then return to 
the turbulence data and explain our measurements of the 
Lagrangian power in terms of Craig's XY-distribution. 



LAGRANGIAN TURBULENCE IN TWO 
DIMENSIONS 

The physics of forced turbulence in 2D differs essen- 
tially from its three-dimensional counterpart. This dif- 
ference can be traced to the presence of a second in- 
viscid invariant in 2D, in addition to the energy. It 
is called enstrophy. The pivotal role of the enstrophy 
was first elucidated in seminal work by Kraichnan[l3j, 
Leith[l3| and Batchelor [l^ (KLB) which has since come 
to be considered as the classical theory of turbulence in 
2D. The essential insight of the KLB theory is that si- 
multaneous energy and enstrophy conservation requires 
the establishment of a second cascade with an indepen- 
dent dissipation mechanism if a stationary state is to be 
reached. In this dual-cascade picture, the enstrophy cas- 
cades from the forcing scale to smaller scales whereas the 
energy simultaneously cascades from the forcing scale to 
larger scales. This latter phenomenon, entirely absent 
in 3D, is known as an inverse cascade. Viscosity dissi- 
pates enstrophy at small scales whereas friction between 
the fluid layer and the substrate on which it moves (Ek- 
man damping) dissipates energy at large scales. Assum- 
ing asymptotically large separation of the forcing and 
dissipations scales for both cascades, and applying stan- 
dard Kolmogorov phenomenology [l6j, the KLB theory 
predicts that the energy spectrum in the direct cascade 
range should scale as fc~^, and in the inverse cascade 
range as k~^/^ where k is the modulus of the wave- vector, 
k. These spectra carry constant fluxes of enstrophy or 
energy through their respective inertial ranges, that is, 
ranges of scales over which forcing and dissipation are 
negligible. This dual cascade theory has been shown to 
be in reasonable agreement with numerical simulations, 
provided that a sufficiently large inertial range is avail- 
able to each cascade[l3- Experimentally, both the direct 
and inverse cascades have been observed While ob- 
serving both cascades simultaneously is rather difficult, 
there is broad agreement that the KLB theory is correct 
asymptotically. See [l^ for a review of experiments. 

Much research has focused on the inverse cascade in 
isolation since it is responsible for the generation of large 
scale coherent motions in two dimensional flow and, de- 
spite the intertwining of the two cascades in the KLB the- 
ory, it is known to persist even when the direct cascade 
range is under-developed[20'|. Numerical experiments 21 
strongly support the k~^/^ scaling for the inverse cascade 



spectrum, laboratory experiments [24] are consistent and 
the phenomenological KLB description of inverse energy 
transfer has been put in a firm quantitative basis [l^. 

We solve the incompressible Navier-Stokes equation 
with Ekman term for the 2D velocity field, v(x, t): 



dv 
'di 



(v • V) V -I- Vp = I'V^v -a^r + t 
V-v 0. 



(1) 



Here, f(x, t) is a forcing term (discussed below) and 
p(x, t) is the pressure field. The parameters v (viscosity) 
and a (Ekman friction coefficient) control the dissipation 
at small and large scales respectively. As a result of using 
physical dissipation terms rather than hyper- viscosity, for 
example, we do not develop a large inertial interval in our 
simulations, which are quite small. The modest inertial 
intervals in our simulations do not pose problems. We 
are interested in studying energy injection into a turbu- 
lent flow. Unlike properties which follow from scaling 
arguments, this is not an asymptotic property and does 
not require R ^ 00. Nevertheless we use serveral dif- 
ferent simulations at different values of R to probe the 
robustness of our results. The simulations are done in a 
doubly-periodic box of the size, L = 27r, using a stan- 
dard pseudo-spectral solver with full dealiasing. For the 
inverse cascade simulations, we used computational do- 
mains of sizes 256^, 512^ and 1024^ in order to investigate 
the effects of varying Reynolds number. For comparison 
purposes, we also performed a small simulation of the di- 
rect cascade regime at a resolution of 256^. Time integra- 
tion was done using a 3rd order Runge-Kutta integrator 
with integrating factors. 

Unlike many numerical simulations of turbulence, the 
forcing term which we used is deterministic. It was de- 
signed to model the electromagnetic forcing which is a 
popular experimental method of driving turbulence in 
thin fluid layers [i^, 24 1. The idea is to place an ar- 



ray of magnets underneath the fluid layer in some par- 
ticular arrangement. One then passes an electric cur- 
rent, which may have some non-trivial time dependence, 
through the fluid layer so that the Lorentz force pro- 
duces a quasi-2D body force on the layer. Suppose that 
the electric current is applied in the x-direction. Denote 
it by l{t) = I{t) x, with x being the unit vector in the 
x-direction. If B(x) = i?(x) z denotes the (vertical com- 
ponent of) the magnetic field generated by the magnet 
array, then the body force exerted on the fluid is 



f(x,i) = g\{t) X B(x) 



(2) 



where is a phenomenological coupling parameter which 
measures the strength of the coupling of the fluid to 
the magnetic field. Having applied the current in the 
x-direction, the force acts purely in the y-direction. In 
this article, we consider both the direct current (DC) case 
in which I{t) is independent of time and the alternating 
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FIG. 1: (Color online) (A) Magnetic field used to force the (256x256) inverse cascade simulation. (B) Typical vorticity snapshot 
from the inverse cascade. (C) Magnetic field used to force the (256x256) direct cascade simulation. (D) Typical vorticity 
snapshot from the direct cascade. 



current (AC) case where I{t) is a sinusoidal. Both are ex- 
perimentally relevant. To produce an inverse cascade, the 
magnetic field is chosen to excite modes at small scales. 
For the direct cascade the magnetic field is chosen to ex- 
cite modes at large scales Fig. [HA) shows the magnetic 
field distribution, B(x), used for forcing the inverse cas- 
cade simulations and Fig. [IjC) shows the corresponding 
field for the direct cascade simulations. The field -B(x) is 
generated by choosing a sum of modes in spectral space 
clustered around some characteristic wavenumber. These 
modes are then assigned random phases and an inverse 
Fourier transform taken to produce a spatially disordered 
forcing field. The characteristic forcing wavenumbers are 
around 32, 64 and 128 for the inverse cascade simulations 
and around 3 for the direct cascade simulation. Note 
that, unlike some turbulence forcing schemes which em- 
ploy a time- varying random forcing field, the disorder in 
our system is quenched. That is to say, once the initial 
random phases are assigned to produce a disordered mag- 
netic field such as those shown in Figs. [UJA) andllJC), 
the magnetic field remains fixed for the duration of the 
simulation. 



Figs. IHB) and [HD) show instantaneous snapshots of 
the vorticity field for the inverse and direct cascades re- 
spectively at 256^ resolution. Both are in the developed 
turbulence regime. The two fields are qualitatively very 
different. The inverse cascade vorticity field contains 
many small incoherent vortices at the scale of the forc- 
ing coexisting with larger scale clusters of like-sign vor- 
tices. The direct cascade vorticity field is dominated by 
a smaller number of more coherent vortices, which are 
again comparable in size to the forcing scale, separated 
by long ribbons. 

The two regimes also differ qualitatively in spectral 
terms. Instantaneous energy spectra are shown in Fig. [2] 
and FiglSl The inverse cascade spectrum is very close 
to the KLB prediction of despite the fact that the 

simulation does not attempt to resolve the direct cascade 
range. It is accepted that the inverse cascade scaling of 
^-5/3 (^Qgg jjot require an extensive direct cascade range 
despite the fact that it is required in the theoretical ar- 
gument [i^l • On the other hand, the direct cascade spec- 
trum is much closer to k~*-^ than k~^. This is again 
in agreement with extensive numerical and experimental 
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(A) Force PDFs (B) Velocity PDFs 
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FIG. 2; (Color online) Snapshots of the energy spectra in the 
inverse cascade regime at resolutions of 256^ (blue *), 512^ 
(red x) and 1024^ ( green o). Inset (A) shows the correspond- 
ing pdfs of the Lagrangian force and associated best fit Gaus- 
sian distributions. Inset (B) shows the corresponding pdfs 
of the Lagrangian velocity and associated best fit Gaussian 
distributions. 
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FIG. 3: (Color online) Snapshot of the energy spectrum in the 
inverse cascade regime at a resolution of 256'^. Inset (A) shows 
the pdf of the Lagrangian force and associated best fit Gaus- 
sian distribution. Inset (B) shows the pdf of the Lagrangian 
velocity and associated best fit Gaussian distribution. 



investigations of the direct cascade where it is typically 
found that the presence of the coherent vortices at large 
scales produces spectra which are steeper than the KLB 
prediction [2^ and the k^^ is only observed when very 
large ranges are considered for both cascades 

When one talks about properties of parcels of fluid 



moving in turbulence, it is natural to adopt a Lagrangian 
perspective. That is, one considers a set of N tracer 
particles, having positions, x^, i — I . . . N, which follow 
the fluid flow passively. The tracer positions satisfy the 
advection equation: 



dt 



= v{x,{t),t) i = l...N 



(3) 



where v(x, t) is the solution of Eq. H]). To study the 
injection of energy into the turbulence from the forcing 
field, we computed the evolution of 100 such Lagrangian 
tracers. The Lagrangian power is defined 



F(t) =v(x,(t),i) •f(x,(i),t). 



(4) 



For our choice of forcing, f is purely in the y-direction 
so that the power is a simple product rather than a 
dot product, P{t) — Vy{xi{t),t)fy{xi{t),t). We collected 
time series of Vy{xi{t),t) and fy{xi{t),t) for each tracer. 
Multiplying these two together gave us a timeseries of 
the Lagrangian power. The velocity in a turbulent flow, 
be it Lagrangian or Eulerian, is a fluctuating quantity. 
Although the forcing fleld in Eq.Q is deterministic, it 
is sampled along a random trajectory followed by the 
particles, so that fy (x^ (t) , t) is also a fluctuating quanti- 
ties. We computed the empirical probability distribution 
functions of both Vy{xi{t),t) and fy{xi{t),t) from our nu- 
merical data. The resulting distributions are shown for 
the inverse cascade simulations in the insets of Fig. [2] and 
for the direct cascade simulations in the insets of Fig. [S] 
Both are close to Gaussian [26]. These results suggest 
that the Lagrangian power can be modeled as the prod- 
uct of two almost Gaussian variables which are presum- 
ably correlated to some degree [sj]. We shall present 
measurements of the power in due course and test this 
assertion. First, however, we need some understanding 
of the probability distribution of products of correlated 
Gaussian variables. This is addressed in the next section. 



PRODUCTS OF NORMAL VARIABLES: 
CRAIG'S XY DISTRIBUTION 



The product of two normally distributed random vari- 
ables was first considered by Craig 27[. Consider two 
random variables, x and y, which follow a normal bi- 
variate distribution with means and fj,y respectively, 
standard deviations ax and ay respectively and correla- 
tion coefficient p. Let Z — xy. In [27| . an expression for 
the moment generating function of Z was derived and 
studied and the distribution function, lixy(Z), was ex- 
pressed as a difference of two integrals. For the purposes 
of the present work, we shall restrict ourselves to the case 



of zero means, jix = l^y = 0, and derive Ilxy{Z) directly 
in a form appropriate for asymptotic analysis. 
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For this case the joint distribution of x and y is [2 
1 ^ 



(5) 

To obtain the distribution of Z we begin with the stan- 
dard construction: 

Ii^y{Z) = E[6ixy-Z)] (6) 
= — / dw e^^™ E [e-*"^*"] . (7) 

Here, the expectation is with respect to the distribution 
Calculating E [e""^™] is relatively simple: 

/oo />oo 
dx / dy n(x, y) e^"^"' 
-OO J — oo 

= K== [ e-i-■4(-)^ (8) 

where x = {x,y) and 



A{w) 



-^{i-P^) 

iW n TT 



Evaluating the integral in Eq. ((8|) is a straightforward ap- 
plication of Gaussian integration (see, for example (29|): 



dx e-5''-'^("')'' = 



(27r) 



AciA{w) 



27r 



where 



A^ 
A- 



1 



0-xCrj^(l + p) 
1 

CT.TCr.y(l - p)' 



(10) 
(11) 

(12) 

(13) 



We are left with the following expression for the product 
distribution: 



nxy(^) = 



VA+F 
2^ 



'^dw 



^y{w ~ iA+){w + iA-) 



(14) 



This integral cannot generally be expressed in terms of 
elementary functions, except in the case of p = 0. In 
this case, A~ = A+ = l/{aj-ay) and Eq.([T4]) becomes 
proportional to one of the integral representations of the 
zeroth order modified Bessel function of the second kind 
M : 



nxy(^) = z^^o(r^). 



nax<Jv 



axCTv 



(15) 
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FIG. 4; (Color online) Craig's XY-distribution for different 
values of p. 
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FIG. 5: (Color online) Integration contours for Eq.(l) 



For general values of p, the integral can be evaluated 
numerically. Some curves are shown for a range of posi- 
tive values of p in FigHl We remark that the pdf is always 
peaked at zero with asymmetric exponential tails. The 
degree of asymmetry increases with increasing p. This 
asymmetry reflects the fact that as the degree of corre- 
lation between x and y increases, they are increasingly 
likely to have the same sign. This ensures that the prob- 
ability of a positive value for the product increases with 
increasing p while the probability of a negative value de- 
creases. For negative values of p, or anti-correlation be- 
tween x and y, the asymmetry is in the opposite sense. 

Although we cannot evaluate Ilxy{z) in general, it is 



6 



possible to calculate the asymptotic behaviour for large 
absolute values of Z for any p. Let us look at the struc- 
ture of the integrand in Eg. p^ in the complex ui-plane. 
The analytic landscape is shown in Fig. O There are 
two singularities on the imaginary axis at lA"*" and — iA~. 
The square root necessitates the introduction of a branch 
cut joining these two singularities which we take to be 
along the imaginary axis at (zA+, zc»] and [— ioo, —ihr\ 
Depending on the sign of we deform the original con- 
tour of integration along the real axis into either the up- 
per or lower complex plane as shown in Fig. [51 Careful 
study of the integrand reveals that it acquires a phase 
difference of tt upon going from one side of the branch 
cut to the other. Considering the case Z > 0, we may 
thus write: 



nxy(^) 



VA+iT 



27r 



' dw 



A+ ^J{w + ik^){w - iA-) 
y/{w + iA+){w " iA-) 



Rescaling w ^ iw and putting the two integrals together 
we get 



n,y(z) = 



VA+F 



dw 



A+ A+)(w + A-) 



(16) 



Changing variables, w = A+ + and performing some 
algebra, we can write this in the form 



n.y(z) = - 

TT 



A+A- 



A+ + A- 



du u 



1 



A-)Z 



(17) 



To finally obtain an asymptotic series, we expand the last 
factor using the Binomial Theorem and integrate term 
by term (recalling the definition of the Gamma function: 
T{x) = / t^'-^e-^dt). We obtain: 



nxy(Z) 



A+A- 



-A+Z 



A+ + A 
1 



(18) 



k=0 



k 



r(fc + -)(A+ + A- 



FoUowing the same approach for Z < 0, taking the de- 
formed contour in the lower half plane, we obtain a very 
similar formula with A+ and A^ interchanged and Z re- 
placed by \Z\. Noting that r{^) — -^/tF, the leading order 
asymptotic behaviour for large absolute values of Z is 
found to be: 



nxy(^) 



A+A- 



-A+z 



7r(A++A-) V2 

A+A- e-^"l^l 
7r(A++A-) ■ 



Z>0 
z <0 



(19) 



The strong cusp at zero is the second striking feature 
of the pdfs shown in Fig. |4l Let us briefly investigate 
the behaviour of the integral in Eq. p^ near Z = 0. A 
fast way to compute the leading behaviour as Z — )■ is 
to make the change of variables w = A"*" -I- m in Eq. , 
then differentiate the resulting expression with respect to 
Z and apply the integrating factor e^^'^, to obtain the 
following differential equation: 



iz)) = - 



Va+iP 



ue 



-Zu 



du 



^Ju{u + A+ + A-) 
(20) 

Putting the integral on the right into Mathematica, we 
find that it can be expressed in terms of a confluent hy- 
pergeometric function of the second kind which has a 
tabulated Taylor expansion for small values of its argu- 
ment fH : 



ue du 



^ [/fi,0,(A+ + A-)Z 



yju{u + A+ + A-) 2Z V2 



2(A+-fA-)Z0F 



0(1). 



Sustituting this back into Eg. ((20)) and integrating term 
by term gives, to leading order in Z: 



Va+a^ 



logZ + 0(l) 



(21) 



Note that this simplistic approach does not provide an 
obvious way of determining the constant which arises 
upon integrating Eq. ([20]) . Thus we do not obtain the 
subleading terms in the expansion near zero. For this, a 
more sophisticated analysis will be required. Neverthe- 
less, this quick calculation is sufficient to demonstrate 
that the product distribution is logarithmically singular 
near zero. Further, we note that, to leading order, the 
singularity is symmetric about zero. 



STATISTICS OF LAGRANGIAN POWER 
(DIRECT CURRENT) 

Let us now return to the Lagrangian power. We have 
seen from the insets of Fig. [2] and Fig. [3] that in both 
the inverse and direct cascade regimes, both the velocity 
and force are close to Gaussian in the Lagrangian frame. 
Some sample time series of the product of the two, the 
Lagrangian power, are shown in Fig. [6l We see that 
both traces exhibit wild fluctuations, both positive and 
negative. The timeseries look far from Gaussian, as one 
might expect from the discussion of the previous section. 

It is clear from Eq. (ll4p that a,j:(7yIl:!^y{axCryZ) is a func- 
tion of p only. Therefore, in plotting the numerical data, 
we use this rescaling together with the values of the stan- 
dard deviations presented in TablelH to collapse all of our 
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FIG. 6: (Color online) Typical timetraces of the Lagrangian 
power in the direct (red solid line) and inverse (blue dashed 
line) cascade regimes normalised by their rms values. 



data to similar curves. The rescaled empirical distribu- 
tion functions of the power for both direct and inverse 
cascades at 256^ resolution are shown in Fig. [T] They 
almost collapse to the same curve and are qualitatively 
similar to Craig's XY-distribution shown in Fig. [4] They 
are strongly peaked at zero with asymmetric exponen- 
tial tails. The positive tail decays more slowly that the 
negative tails which gives the distribution and net posi- 
tive mean value, despite the fact that the most probable 
value is zero. The peak at zero is consistent with the 
logarithmic singularity suggested by Eq. ipT]) but such a 
divergence is too weak to observe unambiguously with 
resolution we were able to achieve for the pdf near zero. 



Simulation 
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IC 512^ 


0.42 


7.3 X 10" 


-5 


3.6 X 10* 


0.063 


0.086 
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IC 1024^ 


0.57 


1.8 X 10" 


'5 


2.0 X 10^ 


0.21 


0.16 


0.027 


DC 256^ 


1.0 


5.8 X 10" 


'i 


1.1 X 10* 


0.34 


0.59 


0.14 



TABLE I: Parameters of the inverse cascade (IC) and di- 
rect cascade (DC) simulations to two significant digits. The 
columns are, from left to right, r.m.s. velocity, Vrms, kine- 
matic viscosity, i/, Reynolds number, R — 2TTVrms/i', vari- 
ance of the Lagrangian force, cr^, variance of the Lagrangian 
velocity, crj, force- velocity correlation coefficient, p. 



FIG. 7: (Color online) Rescaled empirical pdfs of Lagrangian 
power for both the inverse (red -|-) and direct (green x) cas- 
cades at 256^^ resolution. Solid lines indicate the tails pre- 
dicted by Eq. (|19|l for the measured values of the v-f correla- 
tion coefficient, p. 
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FIG. 8; (Color online) Rescaled empirical pdfs of Lagrangian 
power for the inverse cascade at resolutions of 256^ (red 
o), 512^ (green x) and 1024^ (blue *) with corresponding 
Reynolds numbers as noted in Table ID The solid line indi- 
cates the tails predicted by Ea. (|19|l for the most asymmetric 
case (256^) having p = 0.11. 



To make this more quantitative, we should compare the 
empirical distributions with Eq. (fT4|) and Eg. fTS)) . taking 
X to be the force and y to be the velocity. Table U shows 
the measured values of the correlation coefficient, p, be- 
tween the force and velocity and their standard devia- 
tions, (Jf and cr„. Fig. [7] also includes plots of the tails 
predicted from Eq. p9p , for the measured p values shown 



in Table [H In fact, the curves obtained for p = 0.11 
and p = 0.14 are almost impossible to distinguish at the 
level of convergence which we have been able to obtain 
from our data. Notwithstanding the question of why the 
values of p should be so close for the two regimes, both 
show excellent agreement for the tails of the power dis- 
tribution, confirming our qualitative arguments that the 
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pdf should be close to Craig's XY-distribution. 

Next we investigate whether Craig's XY-distribution 
remains a good model of the power pdf as the Reynolds 
number is increased. Fig. [8] shows the rescaled pdfs of 
the power obtained for the three inverse cascade simu- 
lations at resolutions of 256^, 512^ and 1024^. The cor- 
responding Reynolds numbers are shown in Table [H It 
is clear from Fig. [8] that the three pdfs almost rescale to 
the same curve. There is a discernible trend towards de- 
creasing asymmetry as the Reynolds number is increased 
which is reflected in the decreasing value of the corre- 
lation coefficient, p. This is to be expected intuitively. 
As the velocity becomes more turbulent, it becomes less 
correlated with the forcing field. Given that we have es- 
tablished that the average power injected is non-zero due 
to the asymmetry of the tails, one might wonder whether 
the average power also decreases as R increases and the 
tails grow more symmetric. In fact the mean of the XY- 
distribution is UyUfp so a decreasing p can be compen- 
sated for by increasing the variance of the velocity as R 
increases to maintain a constant rate of energy injection. 
Nevertheless, it would be interesting to understand ex- 
actly how p decreases with increasing R. The simulations 
presented here, while sufficient to identify the trend, are 
clearly insufficient to answer this question quantitatively. 

One final point should be made about Fig. [8] which 
clearly illustrates the limitations of this kind of model- 
ing. While the data from the 256^ simulations shown in 
Fig. [7] are almost perfectly fitted by Eq.lUHl), the 512^ 
and 1024^ simulations show a smoothing out of the loga- 
rithmic cusp at zero, although the tails remain quite well 
captured. The XY-distribution only describes the power 
pdf as well as the underlying force and velocity distri- 
butions can be approximated by Gaussian distributions. 
Close inspection of the velocity pdfs for the higher resolu- 
tion inverse cascade simulations indicate that the pdfs of 
velocity become slightly flatter than Gaussian near zero. 
This variation presumably accounts for the smoothing of 
the cusp. One might speculate about the physical mean- 
ing of this deviation but that is not the purpose of the 
present work. Indeed, in the following section, we shall 
see that large deviations from Gaussianity for the forcing 
field can result from simple modifications of the experi- 
mental set up. 



STATISTICS OF LAGRANGIAN POWER 
(ALTERNATING CURRENT) 

Although a Gaussian forcing field of the type observed 
in the direct current simulations is often obeserved ex- 
perimentally, and is often imposed theoretically in order 
to simplify calculations, it is by no means true that the 
forcing must be Gaussian. Thus, while the results de- 
rived above may be applicable to a range of situations, 
we do not claim that they are in any way universal. One 
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FIG. 9: (Color online) Pdfs of the magnetic field (green -I-) 
and velocity (red o) in the Lagrangian frame for the inverse 
cascade regime with AC driving. 



of the simplest ways to produce a non-Gaussian forcing 
is to apply AC current. We shall study this case in this 
section in somewhat less detail than the DC case. The 
principle aim is to compare with the Gaussian case where 
we believe that the simple product distribution discussed 
above provides a good model of turbulent power fluctu- 
ations. For the AC simulations discussed below we take 
the current to be I{t) = cos(27ra;i) with uj = 10. 

For the case of AC driving, the Lagrangian force is 
itself a product of a non-trivial current and a magnetic 
held: 



fy[t) = mB{t). 



(22) 



Let us assume that B{t) remains Gaussian. This assump- 
tion is well supported by the numerical measurements. 
Fig. [5] shows the magnetic fleld and velocity distributions 
for the inverse cascade regime driven by an alternating 
current. Both remain close to Gaussian as in the direct 
current case. So, given that B{t) is Gaussian, how is 
the product f{t) = I{t)B{t) distributed? Of course, to 
deflnitively answer this question we would need to un- 
derstand if and how the magnetic fleld is correlated with 
the current. However the simplest approximation is to 
assume they are uncorrelated, another assumption which 
is supported by our measurements. In this case, we may 
treat I{t) as a random variable generated by uniformly 
sampling the phase. Lot over the interval [0, If we 
proceed with the calculation, this leads after some sim- 
ple manipulations to the following expression for the pdf 
oil: 



MI) = r pe'-'Mw), 
J^oc 27r 



(23) 
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FIG. 10; (Color online) Pdfs ol the Lagrangian force (green 
o) and power (red +) for the inverse cascade regime with AC 
driving. 



where Jq (w) is the zeroth order Bessel function of the first 
kind. This integral has a name - Weber's Discontinuous 
Integral (see [30j, chap. IV), and can be evaluated as 
follows: 



n/(/) 



^/l^ 







if l/l < 1 

if l/l > 1 



(24) 



This distribution makes sense. The current is a cosine 
and cannot take values outside of the range [—1, 1]. Fur- 
thermore, the most probable values are ±1 since this is 
where the cosine function is flatest. Let us now combine 
this distribution with an assumed Gaussian distribution 
for B{t) with mean zero and variance as- Doing some 
lengthy, but elementary, calculations in the spirit of our 
derivation of Ilxy{Z), we obtain the following proposed 
distribution for the Lagrangian force in the case of AC 
driving: 




(25) 



This formula is compared with the measured La- 
grangian force distribution in Fig [TOl We see that the 
alternatic current strongly modifies the distribution of 
the force from the Gaussian distribution measured in the 
direct current case. However, this modification is cor- 
rectly captured by Eq. ips)) . 

Finally, one may ask if the distribution of the power 
may be calculated in this case. Here one encounters a 
problem. We have learned that the correlation coeffi- 
cient between the force and the velocity controls the de- 
gree of asymmetry of the power pdf. In the direct current 
case, where both the force and the velocity are Gaussian, 



the bivariate normal distribution, Eq.([5]), provides a very 
reasonable way to correlate the two. In the alternating 
current case, where the force follows the distribution of 
Eq. ipS]) . it is less obvious how to correlate the two. In 
the absence of a more thorough analysis addressing this 
issue, we give here the result for the uncorrelated case 
only. If the force follows the distribution Eg. ((25)) with 
standard deviation as for the magnetic component and 
the velocity is normally distributed with standard devi- 
ation, (T„, then some lengthy calculations show that the 
uncorrelated product follows the distribution 



H 



(ac) 



iP) = 



1 



dw 



w 



4a| 



(26) 

This distribution is plotted with the measured values of 
as and ay in Fig 1 101 against the empirical distribution of 
the Lagrangian power. As expected, it completely fails to 
capture the asymmetry of the power distribution. Nev- 
ertheless it is sufficiently close to convince us that the 
heuristic arguments presented here are correct and that 
some more careful calculations might be worth doing to 
attempt to incorporate the correlation between v and / 
in a reasonable way. Furthermore, the measured distribu- 
tion for the power looks qualitatively very similar to the 
distributions measured for the direct current case. This is 
despite the fact that the force distributions are so differ- 
ent in the two cases. This suggests that some qualitative 
features of the distribution are more universal than one 
might expect from the very particular calculations per- 
formed above for products of Gaussian variables. These 
questions may be worth further investigation in the fu- 
ture. 



CONCLUSIONS 

We have characterised the injection of energy into a 
two-dimensional flow by measuring the injected power 
in the Lagrangian frame. The measured distribution is 
very sharply peaked at zero with asymmetric exponential 
tails. It can be understood as arising from a simple prod- 
uct of correlated almost-Gaussian variables (the force and 
the velocity). From our analysis, we suggest that the 
power distribution is actually logarithmically singular at 
zero although this is too weak of a divergence to distin- 
guish convincingly in our numerical data. We derived 
an asymptotic expression for the degree of asymmetry of 
the tails of the power distribution which depends only 
on the correlation coefficient of the force and velocity. 
We concluded by comparing the results for a case, alter- 
nating current instead of direct current forcing, for which 
the observed force distribution is far from Gaussian. The 
qualitative features of the power pdf remain unchanged 
although it will be hard to quantify the measurements as 
easily as in the Gaussian case. 
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